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Abstract ^ 

A new hybrid method is presented for the analysis of the scat- 
tering and radiation by conformal antennas and arrays comprised of 
circular or rectangular elements. In addition, calculations for cavity- 
backed spiral antennas are given. The method employs a finite ele- 
ment formulation within the cavity and the boundary integral (ex- 
act boundary condition) for terminating the mesh. By virtue of the 
finite element discretization, the method has no restrictions on the 
geometry and composition of the cavity or its termination. Further- 
more, because of the convolutional nature of the boundary integral 
and the inherent sparseness of the finite element matrix, the storage 
requirement is kept very low at O(n). These unique features of the 
method have already been exploited in other scattering applications 
and have permitted the analysis of large-size structures with remark- 
able efficiency. In this report, we describe the method’s formulation 
and implementation for circular and rectangular patch antennas in 
different superstate and substrate configurations which may also in- 
clude the presence of lumped loads and resistive sheets/cards. Also, 
various modelling approaches are investigated and implemented for 
characterizing a variety of feed structures to permit the computation 
of the input impedance and radiation pattern. Many computational 
examples for rectangular and circular patch configurations are presen- 
ted which demonstrate the method’s versatility, modeling capability 
and accuracy. 
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1 Introduction 


In designing low observable vehicles, the radar cross section (RCS) of con- 
formal microstrip patch or cavity antennas and arrays is of crucial import- 
ance. An example of such an antenna configuration is illustrated in Figure 1 
and in computing its RCS, it is necessary to consider the contributions of 
the radiating elements as well as those of the feed and cavity substructures. 
Conventional analyses based on the space domain and spectral domain in- 
tegral equations [l]— [7] make use of the infinite substrate Green’s function 
and are thus best suited for computing the RCS near resonance. In that 
case, the antenna size is relatively small and structural terminations are not 
as important. Outside resonance and particularly at higher frequencies, the 
moment method solutions lead to large systems which are not practical for 
including contributions due to the feed and substructure of the conformal 
antenna or array. 

In this report, we review a hybrid finite element method for computing 
the scattering of conformal antenna structures while taking into account the 
presence of the feed geometry and the cavity substructure. This method was 
originally presented by Jin and Volakis [8]— [1 1] and has recently been gener- 
alized to allow modeling of non-rectangular patches [12] and more intricate 
feed structures and distributed loads such as resistive sheets [13], often used 
for RCS reduction purposes. An important feature of the method is the use 
of the finite element method to discretize the cavity region occupied by the 
substrate, feed structures and radiating elements. This results in a sparse 
matrix for the solution of the cavity fields and is the primary advantage of 
the method because it leads to 0(N ) matrix system, where N denotes the 
number of edges used in constructing the volume mesh. At the aperture of 
the cavity, the exact boundary integral equation is enforced to terminate the 
mesh (i.e. not an absorbing boundary condition) and thus, in principle, the 
proposed hybrid finite element method is exact. To retain an 0(N ) storage 
requirement for the overall system, it is necessary to have a structured grid at 
the aperture allowing use of the FFT in conjunction with an iterative solver 
[14], [15], [16], 

In the following, we first present the mathematical formulation of the 
method, allowing for the presence of lumped loads, resistive sheets and feed- 
ing cables or microstip lines. The resulting functional is discretized by in- 
troducing tetrahedral or rectangular elements. In particular, edge-based ele- 
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ments are employed which have a number of inherent advantages over tra- 
ditional node-based elements. Specifically, the edge-based elements are di- 
vergenceless within each element, thus, eliminating spurious solutions; they 
maintain tangential field continuity at dielectric interfaces without a need 
to introduce additional constraints; and are suitable for modeling corners, 
where certain field components are singular. 

Using codes based on the presented finite element boundary integral for- 
mulation, we compute the RCS of different configurations which demonstrate 
the accuracy of the method as well as its adaptability in modeling a large 
class of configurations. 


2 Formulation 

Consider the cavity configuration illustrated in Figure 1. The cavity is re- 
cessed in a ground plane as is usually the case with conformal antennas on 
aircraft platforms. The actual antenna or array elements are placed either 
on the aperture or within the cavity, in which case they radiate in a sub- 
strate/superstrate configuration. It is assumed that the elements are fed in 
various ways, including via a stripline or a microstrip line, coaxial cable or 
alternatively, they may be aperture coupled to the feeding lines. Further- 
more, this formulation allows for the cavity material and geometry to be 
arbitrary, and the same is true for the shape of the patches. The patches 
may be further loaded with resistive cards or lumped elements for bandwidth 
or frequency control. For the described cavity-backed antenna configuration, 
we are interested in computing its radar cross section over a wide range of 
frequencies. 

In accordance with the method introduced by Jin and Volakis [11], the 
problem described above can be formulated by separating the computational 
domain into two regions, to be referred to as regions I and II. In region II, 
which encompasses the volume occupied by the cavity, the finite element 
method will be employed for discretizing the fields. The primary reason for 
using the finite element method is its adaptability in modeling a variety of 
cavities and radiating elements whose geometries and material composition 
may vary widely. The fields in region I (exterior region) will be computed 
via the boundary integral method. This amounts to introducing equivalent 
sources on the cavity aperture which are then integrated to obtain the radi- 
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ated fields. The final step of the formulation is to couple the fields in each 
region by enforcing tangential field continuity at the aperture. When the 
final equations are discretized, the resulting system involves a sparse section 
and a full submatrix. The first is associated with the interior fields and, 
thus, involves the bulk of the unknowns. The full submatrix is associated 
with the field components on the boundary and by resorting to a structured 
grid, this submatrix is Toeplitz which is important for retaining the 0(N ) 
storage requirement of the proposed hybrid finite element method. 


2.1 Interior Region Formulation 

The goal with any partial differential equation (PDE) solution method, such 
as the finite element method, is to solve the vector wave equation 

v X f-j-V x E ) - kltr E = -jk 0 Z 0 r fVx|-i-Mj (1) 

where (e r ,/z r ) denote the relative constitutive parameters of the material 
filling the cavity, k 0 is the free space wavenumber and Z 0 is the free space 
intrinsic impedance. In addition, J* and M‘ denote the impressed (if any) 
electric and equivalent magnetic sources, respectively. These are generally 
set to zero for scattering calculations but we shall retain their presence for the 
sake of generality. Instead of solving (1) directly, it is more appropriate for 
numerical purposes to construct its equivalent weak-form. This can be readily 
done via the weighted residual method. Accordingly, (1) is first weighted with 
the function T(r) giving the expression 


< R - T > = //X{ vx fe VxE )- T -**^- T }* 

+ JJJ V T • [,w + V X ^M' j dv = 0 (2) 


in which V denotes the volume occupied by the cavity. It is desirable for 
numerical calculations to make the first volume integral of the weighted re- 
sidual (R,T) symmetric with respect to E and T. To accomplish this we 
recall the identity 




x E) • V x T dv 
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+ jk 0 Z 0 (jj) s T • (H x n) ds (3) 


in which S encloses the volume V , n denotes the outward unit normal to S 
and H = V x ~E/(jk 0 Zo ) is the magnetic field. When (3) is employed into (2), 
we obtain 


<R,T) = 


//Xt (VxE)(VxT) -^ ET )^ 

+ jkoZo ff T [(H 1 + H r + H 5 ) x n] ds 

+ JJI' 


jk 0 Z 0 3 i + V x (— M' 

\f*T 


dv = 0 


( 4 ) 


where the closed surface in (3) has been replaced by the cavity aperture sur- 
face ^ap since E x n vanishes on the interior cavity walls. Consequently, we 
have set H = H' + H r + H S , where H* denotes the incident plane wave excit- 
ation and H r is the reflected field from the ground plane with the aperture 
removed. The remaining field H s is that scattered by the aperture and must 
be related to the electric field on 5 ap either through a local or a global bound- 
ary condition. A local boundary condition, such as an absorbing boundary 
condition, will only provide an approximate relationship between E and H 
on S'ap and is not therefore considered in the context of this formulation. 
In contrast, a global boundary condition such as the Stratton-Chu integral 
equation is an exact relation between E and H on S ap . The treatment of 
the surface integral over 5 ap will be discussed later in more detail. For the 
moment we return to the overall numerical solution of (4) for the field E. 

To generate from (4) an appropriate linear system, it is necessary to first 
discretize the volume V in smaller volume elements. The field in each ele- 
ment can then be approximated by introducing some linear or higher order 
expansion. Subsequently, by using different testing or weighting functions 
T = W j ( j = 1,2,..., AT) a linear set of equations can be generated for the 
solution of the field expansion coefficients. Some typical volume elements 
for discretizing V are illustrated in Figure 2 and any one of these can be 
employed for meshing the cavity volume. Obviously, elements such as rect- 
angular bricks are best suited for rectangular patches, whereas tetrahedrals 
can be used for discretizing cavities which house antenna elements of more 
arbitrary shape. An advantage of the rectangular bricks is their geometrical 
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simplicity, leading to a rather simple geometry preprocessor. In contrast, the 
mesh generator associated with the tetrahedral elements must be far more 
sophisticated. Nevertheless, regardless of the element type used, the fields 
within each element (say the eth element) can be expressed as 

E = f>'W' (5) 

i=l - - - - 

In this, N v denotes the number of edges forming the element (twelve for 
bricks and four for tetrahedrals) and W' are the element basis or expansion 
functions. For our analysis they have been chosen such that the expansion 
coefficients E- represent the value of the field component along the ith edge 
of the eth element when evaluated at the location of the same edge. 

Substituting the expansion (5) into (4) and setting T = W' (Galerkin’s 


testing) gives 

F\v + Fs = Fi + Ff 

(6) 

where 

N c N c N v 

F,v = £ c,V = £ £ 

(7) 


N e N v 

Ff = £*? = £*■ 

e=l e =l 

(8) 


Nm N^ 

Fs = £ F/, F, = £ Q' 

(9) 


€=1 €5=1 

e €-Sap €€-S»p 


with 



II 

JJJv Q-(V x W') • (V x W')- k 2 0 e T W' ■ W') dv 

(10) 

K' = 

-JIL w 'r {' w, '- Vx (H)* 

(11) 

Fs = 

jk ° z °IL s ' j ■ [ rS x ”1 ds 

(12) 

Qi = 

-2jk 0 Z 0 JJ S ' S'- [H* x n] ds 

(13) 
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The number of volume elements comprising V has been denoted by N e 
whereas the number of surface elements comprising 5 ap has been denoted 
by N se . For this application, h = £ and SJ is the same as when the last 
is evaluated on the face of the volume element coinciding with 5 ap . 

If resistive cards, lumped loads or coaxial cables are considered, the left- 
hand side of (6) must be supplemented hy other terms. Specifically, in the 
presence of a resistive card (of resistivity R, occupying the surface S r = 
£5*), we must introduce the function 


F W = M E E E' If i(nxS'I.(»xS-)*= E U 4 ) 

e£Sr 1=1 JJS $ n e€5r i=l 


in which S * denotes the surface that belongs to the eth element, has N sr sides 
and coincides with 5 r . Correspondingly, S* are the W* basis functions when 
evaluated on 5*. When a lumped load or conducting post is placed at the 
zth edge of the eth element, we must then introduce the function 

F 3V = J - k ^° 12 E'W* ■ Wj = E\C\: (15) 

with W' ■ evaluated at the location of the load or conducting post. Usually, 
Cfj = 0 for i 7^ j. In (15) Zi denotes the impedance of the load and / is its 
length, equal to the length of the ith edge. 

Finally, when the antenna is fed with a coaxial cable, we must add to the 
left-hand side of (6) the term (see fig. 3) 


F 4 v 



(N ic 

E E^t 

e€S c * b W=1 




7T p 


ds 


(16) 


As usual, 5' ab denotes the surface over the face of the eth element whose 
face, having N sc sides, coincides with 5 ca b- The basis functions are the 
same as W'j when the last is evaluated on and 7o is the strength of 
the current excitation on the center conductor. Also, p is the radial distance 
measured from the same center conductor and p is the associated unit vector 
along the radial direction. We note that the expression (16) was derived on 
the assumption that the coaxial cable supports the lowest order mode. 


7 



Using the aforementioned matrix elements, (6) can be rewritten in dis- 
crete Formas - 

£ {[^y + m + [cj] + py} {£;} + f s 

e=l 

= EWI + E W E OU 

c=l «- 1 eES c *b 

«€ Sap 

After the element assembly, the above sum results in a sparse matrix which is 
the primary advantage of the method. The sparsity of this matrix is assured 
because the domain of the shape functions is restricted within a single 
element. Obviously, a single edge can be shared by more than one element 
leading to more equations than unknowns which are reduced to the actual 
number of unknowns during the assembly process. 

The evaluation of the matrix elements Ay, By, etc. is straightforward 
once the basis functions W* are specified. The edge-based expansion func- 
tions for rectangular bricks and tetrahedral elements are given in the Appen- 
dices along with the corresponding expressions for the matrix elements. 


2.2 Exterior formulation 

As discussed earlier, to solve the system resulting from (4) or (6), it is neces- 
sary to replace the magnetic field appearing in the surface integral in terms 
of the electric field, i.e. the working variable. An exact relation between the 
tangential electric and magnetic fields at the aperture is 

H = H' + H r ~ JJ S Go(r,r')-[E(r')x£]ds' (18) 

where Go(r, r') is the free space Green’s function 

e ->fc 0 |r— r'| 

47r|r — r'| 

with r and r' representing the observation and integration points, respect- 
ively. Also, I denotes the unit dyad. Substituting (18) into (4) or (12) and 
invoking the divergence theorem twice along with the identity 

V • (G 0 E) = E • VG 0 + G 0 V • E (20) 



G 0 (r, r') 
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it follows that 


Fs = 2kl JJ S ' {z x S*(r)) • [E(r') x z\ G 0 {r, r') ds'J ds (21) 

+ 2 JL Vs ■ ( SjC(r) x £ ) {/X v ' 5 ■ [E(r,) X ^ Go(r ’ r ' } *'} ds 

On the aperture, the expansion (5) becomes 

£M = IXS''(r) (22) 

«=1 

where S* (r) are the basis (r) given in the Appendices when evaluated on 
Si p which is the surface element bordering the e'th volume element. When 
(22) is introduced in (21), we obtain 

N« N s 

Fs = E E ( 23 ) 

e'=l i=l 

with 

= 2 kl jj [f x S*(r)] ■ | If [Sf(r') x i] G„(r,r')<fe'} ds (24) 

+ 2 // s< v s • (SJ(r) x z) | IJ S ', V s ■ [Sf'(r') x i] G„(r, r') ds'} ds 

We observe that the matrix elements Gff will be non-zero even when 5j(r) 
does not overlap with 5*’(r). In contrast, the matrix element A '• given by 
(10) will be zero whenever W' does not overlap with the expansion func- 
tion Wf. Consequently, the matrix [G^ e ] will be full, whereas the matrix 
resulting after the assembly of [A^] will be very sparse and, thus, mostly 
populated by zeros. 

To proceed with the numerical evaluation of the elements Gf. e , it is neces- 
sary to introduce the explicit expressions for the basis functions S'(r). For 
rectangular bricks, these are equal to the basis functions = p- W® (p = x 
or y ) defined in appendix A with the local z variable set to zero. With these 
basis functions, the surface area S* is a rectangle and therefore standard 
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numerical integration routines can be used except when e f = e. In that case, 
G 0 (r, r') is singular and the integrals must then be evaluated with care. 

In accordance with the above, when the eth and the e'th elements are not 
the same or adjacent, then midpoint integration can be used to find 


Gti 


-25* p 5'pG' 0 (r c ,r;) {k 2 0 S'(x c ,y c ) • S ?(x' c ,y') (25) 

+ 25^ p ^pVs • [S](x c , y c ) x z] V s • [sf(x',y') x z] G 0 (r c ,r' c )} 


where and 5® p are the areas of the associated elements, r c = x c x + y c y 
is the midpoint of the eth surface element and likewise r 7 is the midpoint of 
the e'th surface element. When the eth and eTh elements are adjacent or 
coincident, the rectangular surface areas can be subdivided into 3x3 small 
rectangles. On replacing the expansions functions with their midpoint value 
at each of the nine sub-elements we get 


3 3 

G*/ ~ 2 'ff 'ff 5 mnpq { — ^O^j( r ”»n) ■ S, (r p9 ) 

m ,n= 1 p,q = \ 

+ Vs • [S'(r mn ) x z)V s ■ [sf (r p9 ) x £]} (26) 

where r mn and r pq denote the midpoints of each sub-element, and 

*™=JiJJk G ° iT ' r,)ds ‘} ds (27) 


The kernel of g mnpq is singular when r = r', in which case it is appropriate 
to rewrite it as 


9mnpq 



1 

47r|r — r 7 ! 


) 



ds 


+ ff ff — r— r ds' ds 

JJs' mn JJs'f 4?r|r — r'| 

The first of these integrals is now well behaved and can be evaluated nu- 
merically using a 2 x 2 Gaussian integration. The second integral, although 
singular, can be evaluated analytically as described in [17]. 

In the case of tetrahedral elements, the surface basis S'(r) become the 
same as those given by Rao, etc. [18] when defined for a pair of triangles 
which share an edge. The evaluation of the matrix elements G-J for this 
case is more involved and is best accomplished using area coordinates as 
done in [19, 20]. 
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3 Computational and Modeling 
Considerations 

3.1 Mesh generation 

A first step in any numerical analysis is the generation of an appropriate dis- 
crete model of the structure. For our case, this amounts to subdividing the 
cavity into bricks or tetrahedrals, a procedure referred to as mesh generation. 
For rectangular cavities and patches, the mesh generation using rectangular 
bricks is fairly simple. An approach which was employed is to subdivide the 
cavity volume into layers of thickness Az. Each layer is then readily sub- 
divided into A r se x N se rectangular bricks. If the cavity housing the antenna 
structure is nonrectangular, rectangular bricks can still be used, but in that 
case staircasing will be unavoidable. The discretization of the cavity us- 
ing tetrahedrals requires a more sophisticated mesh generation scheme and 
several commercial packages are available for this purpose. We have used 
a package marketed by SDRC IDEAS which can be run on many different 
workstation platforms. 

Regardless of which mesh generation scheme is employed, the final result 
is a file often called the Universal file. As a minimum, this file contains the 
following Tables: 

1. A list of the global volume element Nos. and their material code. 

2. A list of the nodes forming each volume element. 

3. A list of the (x,y,z) coordinates for each node. 

4. A list of global edge Nos. and a list of the pairs of global node Nos. 
forming an edge. 

The last table is not usually contained in the Universal file but can be 
readily constructed from the first three tables. To the above lists we must 
also add information pertaining to the antenna geometry, feeding lines, resis- 
tive cards, lumped loads and about the presence of a coax cable or a probe 
feed. This information is generated by a preprocessor based on information 
provided by the user. For example, in specifying the location of a rectangu- 
lar patch or stripline, the user may enter the z location of the patch along 
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with the (x,y) coordinates for two opposite corners. For non-planar metallic 
surfaces, it is necessary to flag or group the appropriate surface elements 
bordering the metallic surface. This must, of course, be done when gener- 
ating the mesh. A similar procedure is required for identifying the presence 
of resistive cards. The resistivity must be specified for each resistive surface 
element or patch. 

The specification of lumped loads or conducting posts is relatively simple. 
These are assumed to be placed at the edges of the volume elements and 
consequently the mesh must be generated with this in mind. Once the mesh is 
completed, the user specifies the location or edge number and the impedance 
value of each load or post. If the post is perfectly conducting, i.e. a short, it 
is flagged as such and during matrix assembly, the field coefficient associated 
with the edge occupied by the short is set to zero. The center conductor of 
the coaxial cable is defined in this manner. 

Based on the antenna and cavity structural and loading specifications, 
the following Tables are generated: 

5. List of edge Nos. on a perfectly conducting surface (cavity wall, patch, 
stripline, wires, metallic posts, etc.). 

6. List of volume elements along with their global edge Nos. which border 
the aperture and do not coincide with a patch or other metallic surface. 

7. List of volume elements and global edge Nos. bordering a resistive sur- 
face along with the associated resistivity. 

8. List of global edge Nos. along with subglobal edge Nos. for combina- 
tions. 

The preprocessor (see flow chart in fig. 4) which can be run interactively, 
with or without a graphical user interface, depending on which mesh gener- 
ation package is employed. The matrix elements are computed directly from 
the data provided by these tables using l he algorithms outlined in the pre- 
vious section. The excitation vectors (to the right of (17)) are readily filled 
once the user specifies the location and type of source for antenna parameter 
calculations. For scattering calculations {Kf} and {P*} are set to zero with 
{Qf} being the only excitation vector. 
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3.2 Matrix Assembly and System Solution 


Once all input tables are read, the goal is to assemble the matrix system 



[ 0 ] [ 0 ] 1 {E*} _ {bl } 

[o] [G] \ {Ef} - m 


(28) 


where [A] is a very sparse square symmetric matrix, [G] is an N se x N ae 
square full matrix, the column {E„} denotes the field values at the edges 
not bordering 5 ap , and {E„} are the corresponding coefficients for the edges 
on 5 ap . For scattering computations, the column {b^} is set to zero and 
{6^} is equal to { Q e m } as defined in (13). All matrix coefficients in the 
above system are numbered in accordance with the unique global or their 
own subglobal numbering schemes. If the pth global edge lies on a metallic 
surface or coincides with a metallic wire or post, E^ or E? is zero and is, thus, 
not included in the unknown field column. Consequently the computation 
of A pn is skipped. The non-zero matrix elements A mn are computed by 
summing A?,-, B\ j, and D XJ . Note that with this notation, m and n refer 

to global edge numbers, whereas i and j are local indices for the eth element. 
Obviously, an edge may belong to more than one volume element and in that 
case A is constructed by summing the element equations for that edge, ft 
is important to note that only the non-zero elements of the [A] matrix are 
stored, thus, ensuring the O(N) memory demand of the system. 

The [G] matrix is the same as the [G®J ] matrix, once the local indices 
of the last are transformed to global incides. Because this matrix is full, 
it is important to employ uniform gridding to permit the computation of 
the associated matrix-vector products via the FFT. That is, for rectangular 
bricks, all elements of S^p must be chosen to be rectangular and equal. This 
may place some minor restrictions on the separation of the patches but is 
well worth it, because it leads to substantial computational savings. In the 
case of tetrahedrals, we must choose all surface elements to be right triangles 
of equal size. Again, this is not difficult to achieve with the mesh generator, 
but may place some restrictions on the specific shape and separation of the 
elements. 

The [G] matrix is symmetric by its nature, and by employing uniform 
discretization as noted above it is necessary to store only one of the rows of 
each Toeplitz block comprising the matrix. Many solvers can be used for solv- 
ing this system, including LU decomposition or iterative schemes. Because 
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of the sparsity of [A] and the special form of [<?], iterative solvers are best 
suited as they require O(N) storage when the fast Fourier transform is used 
for computing the product [ G]{E*}. Specifically, we employed a biconjugate 
gradient algorithm for symmetric systems. This algorithm is given in Ap- 
pendix C and its storage requirement is t\.25Nv + 10.5Af se complex numbers. 
We have observed rather impressive convergence rates using this algorithm 
(with diagonal preconditioning) for systems having even more than 100,000 
unknowns. For example, for 180,000 unknowns, convergence is achieved in 57 
iterations. The CPU time per iteration on a I1P9000/750 for this system was 
40 sec. Another system of 25,000 unknowns required 66 iterations per angle 
with an average CPU time of 2 sec per iteration. Generally, the sampling 
rate was 12 to 15 elements per linear wavelength. 


4 Numerical Results 


In this section we present radar cross section computations for different 
cavity-backed antenna configurations. Once the fields in the cavity volume 
are computed by solving (28), only the tangential aperture fields need be 
retained for evaluating the antenna RCS (or radiated field). Specifically, the 
scattered field is given by 


HTr) = —jkoYo^——— [[ (^+#)-[E(r0xl]e^ sin<,(l ' cos ^' sin ^ds' (29) 

27rr JJSir, 


where (r, 0, <f>) denote the usual spherical coordinates of the observation point. 
The RCS of the cavity-backed antenna is then given by 


cr pq 


= lim 4?rr 2 

T— ‘•CO 


|H‘(r) • 


(30) 


in which p and q are equal to either 6 or <f> . 

Below we present a number of radar cross section patterns for different 
cavity-backed antenna configurations. The first few computational examples 
are intended for validation purposes, and the remaining expose the reader to 
RCS reduction techniques. Input impedance computation examples are also 
included. 
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Example 1: RCS of a single patch residing in a rectangular cavity 

The geometry of this configuration is shown in Figure 5 along with the com- 
puted backscatter agg RCS pattern. The patch is 1.448" x 1.038" and resides 
on a dielectric substrate having thickness t = 0.057" and relative permittiv- 
ity t T 4.0. The substrate is housed in a 2.89" x 2.10" x 0.057" rectangular 
cavity recessed in a ground plane. The excitation was a ^-polarized plane 
wave at 9.2 GHz in the xz plane, and the analysis was done using rectangu- 
lar bricks. As seen the computed age pattern is in good agreement with the 
measured data [11], We note that for this example the cavity terminations 
play a more important role because the computation frequency is not near 
the resonant frequency of the patch. The measured data were collected by 
placing the cavity-backed antenna configuration on a 5 foot long low cross 
section body [IS]. 

Example 2: RCS of a circular patch near resonance 

The geometry of this patch configuration is shown in Figure 6 along with 
the calculated a$g backscatter RCS as a function of frequency. The diameter 
of the circular patch is 1.3 cm and resides on a substrate 0.0787 cm thick, 
having e T ss 2.33. The excitation is a plane wave incoming at (0, = 60°, 
4>i = 180°), and the computations were done using tetrahedral elements with 
the substrate terminated at a diameter of 1.6148 cm. However, because the 
calculations are near the patch resonance (~ 8 GHz), the substrate termin- 
ations do not contribute significantly to the RCS of the antenna. As seen, 
our computations are in good agreement with the measured data [19] which 
were collected with the patch placed in a cavity having a different periphery. 

Example 3: RCS of a rectangular patch in a rectangular cavity 

This configuration consists of a rectangular patch 0.5 cm x 0.5 cm resid- 
ing on a substrate housed by a rectangular cavity 1 cm x 1 cm x 0.5 cm. 
The substrate’s relative permittivity is 2.33. In Figure 7, we show the agg 
backscatter RCS pattern of this configuration at 3 GHz as a function of 
angle. The incident plane wave is in the xz plane (principal plane) and the 
computations were done using rectangular bricks and tetrahedral elements 
for discretizing the cavity region. As seen, both discretization schemes give 
identical results. 
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Example 4: RCS of a loaded rectangular patch 

A common method for reducing the RCS of a patch antenna is to apply 
resistive loading. Figure 8 shows the relative RCS reduction which can be 
achieved by placing four loads at the edges of the microstrip patch. It is seen 
that as the load’s value is reduced to zero, the RCS as well as the gain of the 
antenna are also reduced. As is well known [5, 20], the dB RCS reduction at 
resonance is twice that of the gain, and this is observed in Figure 8 when the 
feed’s presence is modeled with a matched load (the feed is at x = 2.5 cm, 
y - 1.7 cm). However, the 2 to 1 dB reduction is not precisely observed 
because the resonant frequency of the patch is a function of the load s value. 
The calculation frequency of 1.97 GHz corresponds to the resonant frequency 
with the loads set to 300 Q each. Unfortunately, as shown in Figure 9, the 
predicted RCS reduction is not maintained outside the operating band of the 
antenna. Away from the resonant frequency, the resistive loads play a much 
lesser role since the RCS is dominated by the antenna structure, including 
the cavity terminations. In a later example, we describe a treatment which 
is better suited for RCS reduction outside the first resonance. 

Example 5: RCS of a RAM coated patch 

Another approach for reducing the RCS of the patch is to coat it with radar 
absorbing material (RAM) as shown in Figure 10. In parity with the study 
performed in [5], this lossy coating is only over the extent of the patch, thus 
minimizing radiation losses. It should also be noted that this specific config- 
uration can only be analyzed using a volume formulation and illustrates the 
adaptability of the finite element method. The relative reduction achieved 
with the placement of a specific coating as a function of its thickness t is 
illustrated in Figure 10. These calculations were performed at the reson- 
ant frequency of the patch which varied from 1.944 GHz at t = 0 cm to 
1.893 GHz at t = 0.3 cm. As a result, the typical relation between gain 
and RCS reduction is observed with much more fidelity. Unfortunately, as 
seen from Figure 11, the observed RCS reduction is again maintained only 
near resonance. At higher frequencies the thickness of the coating causes the 
overall RCS of the structure to be increased. This occurs primarily because 
of its larger electrical thickness. 
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Example 6: RCS of a patch with distributed resistive loading 

With the motivation of reducing the patch RCS at frequencies outside its 
operational band without substantial compromise in gain, we consider the 
placement of a narrow resistive ribbon (R-skirt) around the periphery of 
the patch as illustrated in Figure 12. Because the ribbon’s electrical width 
increases with increasing frequency, it is expected to have a rather noticeable 
effect on the RCS of the patch at higher frequencies. Thus, greater RCS 
reduction can be attained at higher frequencies without compromising the 
antenna’s gain performance. If a certain level of RCS reduction is desired 
at resonance, lumped loads can still be used in conjunction with the R-skirt 
loading. 

The a$e RCS of the patch as a function of frequency depicted in Figure 12 
is given in Figure 13 before and after loading with the R-skirt. The given 
backscatter RCS data correspond to a plane wave excitation incident in the 
xz plane and at an angle of 0 = 70° with respect to the z axis. However, 
as shown in Figure 13, the RCS reduction is fairly uniform over all 6 angles. 
As seen, the resistive skirt loading reduces the RCS in and out of the band 
as opposed to the narrowband performance of the lumped loads. Not sur- 
prisingly, the placement of the R-skirt is seen to slightly change the resonant 
frequency of the patch from 1.9522 GHz down to 1.744 GHz. From the inset 
in Figure 12, the RCS reduction at resonance is about 9 dB for all angles 
of incidences and thus the corresponding gain reduction is approximately 
4.5 dB. 

Example 7: RCS of a cavity-backed spiral 

Cavity-backed spirals are popular broadband antennas whose RCS is cer- 
tainly of concern because of the complex nature of the antenna. Figure 14 
shows the top and cross-sectional view of such an antenna. The spiral is rect- 
angular in shape for modeling convenience and is housed in a metallic cavity 
2.8125 cm x 2.8125 cm with and without absorber at the base of the cav- 
ity. From the given dimensions, the lower limit of the antenna’s operational 
frequency is approximately 3 GHz. We shall therefore consider the RCS of 
the antenna at 4.5 GHz and our goal is to study the effect of the absorber 
at the cavity base and that of a resistive sheet placed uniformly over the 
cavity aperture. The chosen absorber consists of three uniform layers and 
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was designed to be similar to a commercial absorber giving 18 dB reduction 
at 4.5 GHz. 

The principal plane agg backscatter RCS of the spiral antenna using dif- 
ferent treatments is shown in Figure 15 as a function of angle. The cor- 
responding radiation patterns are given in Figure 16. As seen the absorber 
alone reduces the RCS of the antenna by approximately 6 dB and this was 
to be expected (the gain reduction was 3 dB). Clearly, the R card has a more 
dramatic effect on the RCS of the spiral antenna. Some RCS reduction is still 
observed even when the resistivity of the sheet is 10Zo or 5Zo- However, the 
gain of the antenna (not shown) drops substantially more and in particular 
nearly a 10 dB gain reduction was observed with R = 5 Z 0 . For lower val- 
ues of the resistivity the gain reduction is so substantial (as is also the case 
with the RCS of the spiral) making such a treatment impractical. Perhaps 
a better treatment (which reduces the RCS without substantial compromise 
in gain) would be to apply a lossy dielectric coating over the spiral/cavity 
aperture. This coating would be equivalent to a resistive sheet having a res- 
istivity which is a function of frequency. By controlling the thickness and 
loss tangent of the coating, a frequency response may be achieved which is 
appropriate for a specific application. 

Example 8: Input impedance of a rectangular patch 

The specific antenna geometry is shown in Figure 17 consisting of a patch 
3.4 cm x 5.0 cm residing at the aperture of a 7.5 cm x 5.1 cm cavity. The 
substrate is of thickness t = 0.08779 and has a relative permittivity e r = 2.17 
with a loss tangent of 0.0015. Figure 17 shows the complex input impedance 
as a function of frequency with the probe feed placed at (1.22 cm, 0.85 cm). 
These calculations were done in the presence of a 50 f) resistor (lumped load) 
placed between the patch and substrate at (—2.2 cm, —1.5 cm). Our calcula- 
tions are seen to be in good agreement with the measured data [11] whereas 
corresponding calculations using traditional methods are not as accurate. We 
remark that the 50 0 load (50 Q resistor) serves to reduce the RCS of the 
patch at resonance. However, it also reduces the antenna gain. It has been 
observed though that at resonance the dB reduction in RCS is twice that 
of the gain (i.e. if the gain is reduced by 5 dB, the RCS will be reduced by 
approximately 10 dB). Unfortunately, as noted in a previous example the 
resistors do not have any appreciable effect on the RCS when the patch is 
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illuminated with a frequency out of its resonant frequency. 

Example 9: Input impedance of a circular patch 

The specific antenna configuration is a circular patch having a radius of 
1.3 cm and residing at the center of a substrate filled cavity of radius 2.115 cm. 
The cavity depth is t — 0.406 cm (~A/15) and the substrate’s relative dielec- 
tric constant is c r & 2.9. We computed the input impedance of this patch 
placed on the surface of the cavity for a probe feed located 0.8 cm from the 
center of the patch. The computed and measured impedance are depicted 
in a Smith chart shown in Figure 18 by varying the frequency from 3 GHz 
to 3.8 GHz. The agreement between our calculations and reference data is 
indeed quite satisfactory. 


5 Conclusion 

In this report we reviewed a hybrid finite element-boundary integral method 
for computing the RCS of planar antennas. The pertinent finite element and 
boundary integral formulations were presented along with the method’s at- 
tributes. It was demonstrated that by virtue of the finite element method, the 
formulation is adaptable and rather efficient in modeling differing geometries 
and material compositions, including lumped and distributed loads. 

General-purpose computed codes were written on the basis of the presen- 
ted formulation. These were used to generate the presented RCS computa- 
tions for circular and rectangular patch configurations, spiral antennas and 
other non-rectangular configurations. Among the presented examples some 
served for validating the method whereas others exposed the reader to RCS 
reduction schemes. 
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Figure 2. Geometries for various finite element shapes. 
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Figure 4. Computer code flow chart. 
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Figure 6. Comparison of measured and calculated backscatter RCS data for the shown 
circular patch as a function of frequency (#, = 60°, <f>, = 180°). 



= 0.0787 cm 


29 




RCS in dBsm 



Load Resistance in Ohms 

Figure 8. Relative normal incidence RCS/Gain Reduction in dB of the shown rectan- 
gular patch at / = 1.97 GHz as a function of the four lumped loads placed 
at the edges of the patch. The cavity is 0.17558 cm deep and is filled by a 
substrate having a relative dielectric constant of 0.17. 
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Figure 9. RCS of the patch and cavity shown in Figure 2 as a function of frequency. 
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Figure 10. Relative normal incidence RCS/Gain of the unloaded antenna configuration 
in Figure 2 where the patch is coated with RAM having e r = 10.2 — j 3.8 
and /i r = 2.12 - j 1.5. Calculations are referred to the resonant frequency of 
the coated patch. 
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Lower Left corner of patch at (x,y)=(2.1875cm ,2.96875cm) 

Patch dimensions 5.00cm x 3.4375cm 

I F 1 Feed at (x,y) = (3.4375cm,3.75cm) 


Figure 11. Geometry of a patch in a cavity. The patch is loaded with a square resistive 
ribbon having a resistivity of R = 0.5 Z 0 . 
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Figure 13. o$ 6 backscatter RCS of the patch in Figure 11 as a function of angle with 
and without resistive loading. 
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Frequency 4.5 GHz 


Coordinates of cross (X at lower right comer of conducting strip): (1.406 cm, 1.406 cm) 



Spiral with absorber 

/ 


cm 

18 dB Absorber 
at 4.5 GHz 


Figure 14. Spiral configuration with small cavity aperture. 
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Spiral+Absorber 

Spiral+Absorber+Rcard over the aperture (R=5.0) 
Spiral+Absorber+Rcard over the aperture (R=l .0) 


0 in degrees 

Figure 15. RCS of the spiral configuration in Figure 14 


Spiral Zin=109-j340 

Spiral+Absorber Zin=78-j437 

Spiral+Absorber+Rcard (R=5.0) Zin=61-j389 
Spiral+Absorber+Rcard (R=1.0) Zin=31-j395 


0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 90.00 

0 in degrees 

Figure 16. Radiation power pattern of the spiral antenna shown in Figure 14. 
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Figure 17. Comparison of computed and measured input impedance for a loaded mi- 
crostrip patch antenna, (a) Real part, (b) Reactive part. 
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1.0 


Figure 18. Comparison of input impedance computations for a circular patch. The 
patch radius is 1.3 cm and resides in a cavity of radius 2.115 cm and 0.406 cm 
deep filled with a dielectric of c r « 2.9. The feed is 0.8 cm from the center 
of the patch and the frequency is swept from 3 GHz to 3.8 GHz. 
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Appendix A: 

Matrix elements for rectangular bricks 


Referring to Figure Al, it is convenient to first introduce the definitions 

{W*} = {iN^yN^N^Y (Al) 

where j = 1,2, 3, 4. Thus, for example, E\ = (j> x 2, E~ — <j) y 3, = xN x2 , 

W? = yN y3 , etc. The field expansions can then be written as 


Er = EN e xj (y,z)<i>l } -, E y = ^>2 Nyj(z, x) E z = ^2 N*j[x, y) <j> 

j=\ j= 1 i=i 


(A2) 


where 


7 V e — ( fe ~v')( c — g ') • yy e — nlhziH- 7v e — itzklhl - jV e — ^ - • 

J xl — be ' 1 x2 — be > 1 x3 — 6c ’ J x4 — 6c ’ 

yye _ (c-z')(q-r') , yye _ ^(a-x') ■ iye _ (c-z')x' . yye _ *V. ( A31 

JV »1 - ca > iV »2 — ca > iV y3 ~ ca > iV y4 ~ ca » 

yye (a— x')(6— 1 /') . yy e *'(*-»') . Are _ («-*>' ■ Are £^/ 

2 I at J ^2 a6 7 ^3 at * z 4 at 


In these, (x',y\z') denote the local coordinates specifying a point within 
the eth element and from an examination of the expansion functions we 
observe that <f> e xl represents an average of the field component E x along the 
edge segment (1,2). Likewise, </>* 2 is associated with the E x component 
along the edge (3,4) and so on. The elements of the 12 x 12 [A e ] matrix 
defined in (10) can now be readily evaluated in a straightforward manner 
since N PJ ( p = x,y,z) are simple linear functions. To do so, for notational 
and programming convenience, we shall denote these elements as A { p j),( q k), 
where p = x,y,z\ q = x,y,z ; j = 1,2, 3, 4; and k = 1,2, 3, 4. We next define 
the matrices 


[A'o] 


' 4221 ' 
a be 2 4 12 
36 2 1 4 2 

12 2 4 
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' 2 

-2 

1 

-1 ’ 

[*.] = 

-2 

2 

-1 

1 

1 

-1 

2 

-2 


. -1 

1 

-2 

2 . 


' 2 

1 

-2 

-1 ' 

[A' 2 ] = 

1 

2 

-1 

-2 

-2 

-1 

2 

1 


. -1 

-2 

1 

2 . 


’ 2 

1 

-2 

-1 " 


-2 

-1 

2 

1 

[Aal = 

1 

2 

-1 

-2 


-1 

* ) 

L> 

1 

2 . 


where a, b and c are the side lengths of the brick as defined in Figure Al. 
On using this notation we have 


[A] 




[Axj } yk\ 


\Ayj,xk, 




[ A z j'Xk] 


[A z j % zk. 


[v4 r j >x fc] 

[Ayj'Xk] [Ayj,yk] \Ayj,zk\ 

[A z j t xA:] [A z j,z/:] _ 

= _ 6^ |A3 ' T 
_ 6^ |A ' 3|T ' = ~6^7 |/ ' 31 

-^[A's], 1*4+] = -£l Ks) T 

i{^ lA ' 1 + S |A2, - fc “ £ ' MA ' 01 } 


If a resistive sheet of resistivity R borders a face of the eth brick, then Bfj 
will be non-zero. Denoting the field components at the edges of the subject 
rectangular face by E * n E* 2 , E e ]X , and Ej 2 , it follows from (14) that 

jk 0 Z 0 


~ Bi2,i2 


= Bh,H 


= B' 2 J2 = 


3i? 


didj 
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^iU2 — il ” — ^2,il — “*“j 


where dj and are the corresponding side lengths of the rectangle and R 
is the (unnormalized) resistivity of the sheet. All other elements of the [B] 
matrix are zero. 
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Appendix B: 

Evaluation of matrix elements for tetrahedrals 

Referring to fig. B1 and the associated table, the fields in the eth tetrahedron 
are expanded as 

E = £ E‘ w; 


where the basis functions W- are given by 

W e 7 _,(r) = 

f f 7 _, + g 7 _,- x r r <E V £ 

\ 0 outside element 

f 7 -, = 

r, x r, 2 r,j , r, 2 : position vectors of vertices 

^ e ii and ii (see Table) 


bib 

g7 -i = 

6V e 6 ‘ 

P . 

(l't 2 ) 

“l 

6, 

bi = 

|r, 2 - r,,| = length of the zth edge (see Table) 

% = 

element’s volume 

We note that 

V • W? = 0 VxW' = 2g, 


indicating that Wf are divergenceless. Furthermore, 


WJ ( r J ) ■ ej — ; 8, : — | q i ^ j 

where r J has its tip on the jth edge of the tetrahedron. This last property 
ensures that the coefficients Ef = E • e, represent the average field value at 
the ith edge of the tetrahedron. 

Using the above basis functions, we now proceed with the derivation of 
the matrix elements A e tJ . We have 

m — (V x W') ■ (V x W') = f g, • g> V, 

JJJV c l*r 
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Also 


£r Jr. W ' ' W ' il ' = £r III. (* f ' ’ f J> + < r ' D) + ( g . x r ) ' ( Sj x r )} iv 

= e r (/i + h + h) 

where 

D = (f, x gj ) + (fj x g t ) 

and 

/, = /// r, f 

'* = //l r-D* 

73 = JJJv ^ S ‘ X ^ ' ( Sj X ^ ^ 

Since f is a constant vector, /! reduces to 

/l = fi * fiK 

To evaluate /2 we first set 

4 

X = *52 LiXi 

2=1 
4 

y = J2 L >y> 

2 = 1 

4 

2 - 1,2, 
i=l 

where j/,, 2 ; (i = 1, . . . , 4) denote the y, z ) coordinates of the tetrahed- 
ron’s vertices and Li are the simplex coordinates or shape functions for the 
same element. That is, Li is the normalized area of the tetrahedron formed 
by its three corners other than the zth, and the point (x,y, z) located within 
the tetrahedron. Using the standard formula for volume integration within 
a tetrahedral element and simplifying, we have 

h = y D x x t + D y ^2 Vi + D zJ2 Zi 

4 L 1=1 i=i t=i . 
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where D m is the mth component of D. The evaluation of I 3 can be simplified 
by the use of basic vector identities. We have, 

h = g, • gj / |if dv - l (g, • r)(g i • r) dv 
JV e jv* 

= (dtydjy 9iz9jz) / x dv 

J v c 

+ (gix9jx + 9iz9jz) J v y 2 dv 

“I" “f" 9*y9jy) I z dv 

J v e 

~~ {9ix9jy 9jx9iy) I zy dv 
J v c 

(9\ x9jz “I" 9jx9iz ) I zx dv 
J V e 

- (9iy9jz +9jy9i-) I yzdv 
J V c 


where gi m represents the mth component of the vector g,. Each of the above 
integrals can be easily evaluated analytically and the result can be expressed 
in the general form 


J , 


aia r 


20 


mi + £ On ^ ^ O m { 




i=l 


*=1 


for /, m = 1, ... ,3. The parameters a; or a m can represent any of the recti- 
linear variables x, y, z. 
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Appendix C: 

Biconjugate gradient algorithm for solving the 
symmetric system Ax — b 

Initialize the residual and search vectors with an initial guess xo: 
Po = r 0 = b - Ax 0 


Iterate for k = 0, 1 , 2, . . . : 


a k 

x k + 1 
r k+ 1 

Pk 

Pk + 1 

Terminate when 


(. Pk , Ap k ) 

%k + OtkPk 

r k - a k Ap k 

(i'k+i,r k+ i) 

(rk, r k ) 
r k+l + fikPk 


err = 


!b±i!l 

il&li 


< tolerance. 


In the algorithm, (x, y) = x T y, where T denotes the transpose of the column. 
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Figure Bl. Edge and vertex numbering for the tetrahedron. 



